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We report a measurement of the longitudinal double-spin asymmetry All and the differential cross 
section for inclusive midrapidity jet production in polarized proton collisions at y^s = 200 GeV. The 
cross section data cover transverse momenta 5 < pt < 50GeV/c and agree with next-to-leading 
order perturbative QCD evaluations. The All data cover 5 < pt < 17GeV/c and disfavor at 98% 
C.L. maximal positive gluon polarization in the polarized nucleon. 
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Deep-inelastic scattering (DIS) experiments with po- 
larized leptons and polarized nucleons have found that 
the spins of quarks and antiquarks account for only about 
25% of the nucleon spin The gluon helicity distribu- 
tion and orbital angular momenta are thus essential to 
the understanding of the nucleon spin. Analyses of the 
scale dependence of the inclusive nucleon spin structure 
function 0] and recent semi-inclusive DIS data 0] have 
coarsely constrained the possible gluon spin contribution. 
Complementary measurements with strongly interacting 
probes [H, H| give sensitivity to gluons predominantly via 
quark-gluon and gluon-gluon scattering contributions Q . 

In this Letter we report the first measurement of All 
for inclusive jet production in polarized proton collisions, 

ALL + (1) 

where <j ++ and cr H are the inclusive jet cross sections 
when the two colliding proton beams have equal and op- 
posite helicities, respectively. In addition we report the 
inclusive jet differential cross section. 

In pQCD the (un-)polarized jet cross section involves 
a convolution of (un-)polarized quark and gluon distri- 
bution functions and the (un-)polarized hard partonic 
scattering cross section We compare next-to- 

lcading order (NLO) pQCD calculations with the mea- 
sured cross section to test their applicability and to sup- 
port their use in constraining the polarized gluon distri- 
bution through measurement of All- Our data on All 
are sensitive to gluon polarization for momentum frac- 
tions 0.03 < x < 0.3. 

The data were collected at the Brookhaven Relativistic 
Heavy Ion Collider (RHIC) with the Solenoidal Tracker 
at RHIC (STAR) [U in the years 2003 and 2004 using pro- 
ton beams of 100 GeV energy. Typical luminosities were 
2-5 x 10 30 cm~ 2 s _1 . Spin rotator magnets upstream and 
downstream of the STAR interaction region (IR) rotated 
the proton beam spins from and to the stable vertical 
direction in RHIC to provide collisions with longitudinal 
polarizations 0]. The helicities alternated for successive 
bunches of one beam and for successive pairs of bunches 
of the other beam. Thus STAR recorded collisions with 
all beam helicity combinations in rapid succession. 

The polarization of each beam was measured for 
each beam fill with RHIC Coulomb- Nuclear Interference 
(CNI) proton-carbon polarimeters 0, which were cali- 
brated in situ using a polarized atomic hydrogen gas- 
jet target jT^]. Proton beam polarizations were 30%- 
45%. Nonlongitudinal beam polarization components at 
the STAR IR were measured continuously with local po- 
larimeters [IH and were no larger than 9% (absolute). 

The STAR detector subsystems [|| of principal inter- 
est here are the time projection chamber (TPC), the bar- 
rel electromagnetic calorimeter (BEMC), and the beam- 
beam counters (BBC). The TPC tracks charged particles 



in a 0.5 T solenoid magnetic field for all azimuthal an- 
gles ((/>) and pseudorapidities < 1.3. The BEMC is a 
lead-scintillator sampling calorimeter that limited the ac- 
ceptance in 2003 and 2004, covering all 4> and < 77 < 1 
with respect to the TPC center. The BBCs are composed 
of segmented scintillator annuli that span 3.3 < \rj\ < 5.0 
and measure the proton beam luminosity and transverse 
polarization components. 

Proton collision events were identified by coincident 
signals from at least one BBC segment on either side 
of the IR. The cross section for the BBC coincidence re- 
quirement is 26.1±2.0mb, which is 87% of the non-singly 
diffractive pp cross section [12]- The jet data were col- 
lected with a highly prescaled minimum bias (MB) trig- 
ger, requiring a proton collision event, and a high tower 
(HT) calorimetric trigger condition. It required, in ad- 
dition, a signal from at least one BEMC tower of size 
A77 x A(j) — 0.05 x 0.05 above a transverse energy (£r) 
threshold of 2.2 GeV in 2003 (2.2-3.4 GeV at rj = 0-1 in 
2004). In total 2.1 x 10 6 MB and 3.0 x 10 6 HT events 
were analyzed. The integrated luminosity J Cdt amounts 
to 0.18 (0.12) pb _1 for the analyzed 2003 (2004) data. 

Jets were reconstructed using a midpoint-cone algo- 
rithm [13j that clusters reconstructed TPC tracks and 
BEMC energy deposits within a cone in 77 and (j> start- 
ing from energy seeds of at least 0.5 GeV. A cone radius 
fcone = 0.4 was chosen because of the limited BEMC 77 
acceptance. Particle tracks with pt > 0.2 GeV/c were 
considered if they originated from the primary interac- 
tion vertex, which was required to be on the beam axis 
and within 60 cm from the TPC center to ensure uniform 
tracking efficiency. Calorimeter towers were considered 
if their E^ exceeded 0.2 GeV after correction for charged 
hadron contributions determined from TPC tracking. A 
charged pion (photon) mass was assumed for tracks (tow- 
ers) in relating energy and momentum. Jets were re- 
quired to have a reconstructed jet pt > 5 GeV/c and, as a 
tradeoff between acceptance and effects from acceptance 
edges, a reconstructed jet axis intersecting the BEMC at 
nominal n between 0.2 and 0.8. A minimum TPC contri- 
bution to the jet energy, E T pc/E tot > 0.2 (0.1) in 2003 
(2004), was used to suppress apparent jets from beam 
background. The jet resolution was determined to be 
^25% from the momentum balance of dijet events and 
from simulation, and motivated the choice of binning. 

Figure[T]compares the measured and simulated jet pro- 
file \I/(Ar, r CO ne,pr), defined as the average fraction of 
jet Et inside a coaxial inner cone of radius Ar < r conc , 
for the MB and HT data separately. The reconstruc- 
tion software imposed the same trigger requirements as 
in the data. More than 95% of the total jet energy is 
expected to be contained within r conc = 0.4. The HT 
trigger, providing increased selectivity for jets, causes a 
Pt dependent bias toward jets with hard fragments that 
produce an electromagnetic shower. The ^(Ar, rconc, Pt) 
distributions are well reproduced by PYTHlA-based (v 
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FIG. 1: Jet profile *I>(Ar, r cono , px) versus inner cone size Ar 
at r cono = 0.4 for MB (open squares) and HT (filled circles) 
data compared with STAR Monte Carlo simulation in two jet 
p T bins (a) 5.0 < p T < 6.2 and (b) 14.1 < p T < 17.3 GeV/c. 
In (b) the MB jet yield was too small to measure. 

6.205 [HI 'CDF Tune A' settings 0) Monte Carlo simu- 
lations passed through GEANT-based [3| STAR detector 
simulation. The simulations are used in determining the 
cross section and to assess effects of the trigger bias on 
All. In the cross section analysis of HT data an Et 
threshold of 3.5 GeV was imposed on the BEMC trigger 
tower to ensure a uniform trigger efficiency. 

The differential inclusive cross sections were deter- 
mined separately for the MB and HT data according to 



1 d 2 a 



1 N- 



jets 



(2) 



where Nj cts denotes the number of jets observed within a 
pscudorapidity interval Ar; and a transverse momentum 
interval Apt at a mean jet pt- The correction factors 
c(pt) were determined from simulation, and are defined 
as the ratio of the number of jets reconstructed within 
a given pt interval in the simulated data to those gen- 
erated in the PYTHIA final-state particle record. They 
change monotonically for HT events from 0.02 at pt = 
8.3 GeV/c to 0.79 at pt = 43 GeV/c, whereas they are a 
constant 0.69 for MB events with pt < 12.6 GeV/c. Con- 
sistent values were obtained with the herwig [17| gener- 
ator. Typically 35%-40% of the jets generated in a given 
Pt interval were reconstructed in the same interval. Re- 
constructed pt was found to be on average ~20% larger 
than generated p T in each reconstructed p T interval, and 
the difference is taken into account via c(pr)- 

The MB differential cross sections extracted from 1.4 x 
10 3 jets collected in 2003 and 1.1 x 10 3 in 2004 are in 
good agreement (x 2 /ndf = 0.8). A 20% systematic off- 
set for all pt was found between the HT differential cross 
sections extracted from 43 x 10 3 and 42 x 10 3 jets col- 
lected in 2003 and 2004. We ascribe this difference to 
5% uncertainty (included in the systematic errors below) 
in the year-to-year absolute scale of the BEMC calibra- 
tion, which was changed by a factor of ~ 2 between the 
two years, and to uncertainty in the modeling of tem- 
porary BEMC hardware malfunctions. The calibration 
used 20 x 10 6 d+Au collision events in 2003 and 50 x 10 6 
Au+Au events in 2004. The absolute energy scale was 
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FIG. 2: (a) Inclusive differential cross section for p+p — * jet + 
X at ^J~s — 200 GeV versus jet px for a jet cone radius of 0.4. 
The symbols show MB (open squares) and HT (filled circles) 
data from the years 2003 and 2004 combined. The horizontal 
bars indicate the ranges of the j>x intervals. The curve shows a 
NLO calculation [g]. (b) Comparison of theory and data. The 
band indicates the experimental systematic uncertainty. The 
upper (lower) dashed line indicates the relative change of the 
NLO calculation when it is evaluated at pi = px/2 (fj. = 2px). 



set by matching BEMC energy to TPC track momentum 
for well-contained showers from 1.5 < p < 8GeV/c elec- 
trons identified in the TPC. Uncertainties arise in the 
electron selection, from residual hadronic contamination, 
and from the limited d+Au statistics. 

Figure Ufa) shows the arithmetic average of the 2003 
and 2004 MB and HT cross sections versus jet pr- The 
MB and HT data are in good agreement for overlap- 
ping jet pt (x 2 /ndf = 1.0), despite the very different 
c(pt)- The curve shows the NLO pQCD cross section 
of Ref. 0] evaluated at equal factorization and renormal- 
ization scales, fj, = (If = — pr, using the CTEQ6M 
parton distributions [l^ . Figure [2Kb) compares data and 
theory, showing satisfactory agreement over 7 orders of 
magnitude. The theoretical cross section changes by less 
than 23% if /i is varied by a factor of two and increases by 
1% (13%) at pt of 10 (40) GeV/c if the CTEQ6.1M dis- 
tributions are used. The experimental systematic uncer- 
tainty amounts to 8% in the normalization with the BBC 
and 48% in the measured yield, consisting of 5% due to 
residual beam background, 13% on c(p-r), and 46% from 
a 9% uncertainty on the jet energy scale. The BEMC cal- 
ibration and undetected neutral particles dominate in the 
latter. No corrections were made for the nonperturbative 
redistribution of energy into and out of the jet by the 
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underlying event and out-of-cone hadronization. We es- 
timate that such corrections would increase the measured 
differential cross section by ~25% for p T > 10 GeV/c. 

The asymmetry All was extracted for 5 < pt < 
17GeV/c from a HT data sample of about 110 x 10 3 
jets in 2003 and 210 x 10 3 in 2004. The sample size is 
larger than in the cross section analysis, since no BEMC 
energy threshold was required. The jet yields N were 
sorted by equal (++) and opposite (H — ) beam helicity 
configurations. The asymmetry was extracted as: 



All = 



E(PiP2)(N++ - RN+-) 
J2(P 1 P 2 ) 2 (N++ + RN+-) 



(3) 



where P\,i are the measured proton beam polarizations, 
R ~ 1.1 is the ratio of measured luminosities for equal 
and opposite proton beam helicities, and parity violating 
differences < O(10" 4 ) in the cross sections for different 
beam helicities are not considered. The sums are per- 
formed over runs typically lasting 20 minutes. 

The results for All from 2003 and 2004 data are in 
good agreement (x 2 /ndf = 0.3). Figure [3] shows the 
combined All versus jet pr> together with the statis- 
tical (bars) and systematic (bands) uncertainties. 

A 25% combined scale uncertainty arises from the CNI 
beam polarization measurement (22% in 2003 and an un- 
corrected 16% in 2004) and from the CNI absolute cali- 
bration (18% common to both years). 

The uncertainty in R was estimated to be 0.003 using 
narrow and wide timing requirements for the BBC coinci- 
dence. It takes into account differences in sampling of the 
longitudinal vertex distribution in the jet analysis and in 
the relative luminosity measurement, and corresponds to 
0.009 uncertainty in All- An independent measurement 
with the zero degree calorimeters (ZDC) [8j gave consis- 
tent results to within statistical uncertainties. No double 
helicity asymmetry of the BBC measurement relative to 
the ZDC measurement was observed. 

Residual nonlongitudinal proton beam polarization 
at the STAR IR could contaminate the All measure- 
ment through an azimuthally uniform two-spin asymme- 
try [ijj]. A limit of 0.010 on such contamination was set 
from local polarimetry data and from two-spin asymme- 
try measurements with vertically polarized beams. 

Beam background occasionally caused BEMC signals 
not associated with collisions at the IR. Its effect on the 
jet yields was reduced with the aforementioned selection 
on Etpc /-Etot • Residual yields were estimated to be no 
larger than 8% (5%) in the 2003 (2004) data from the 
variation of jet spectra with beam-background conditions 
monitored with the BBCs when filled and empty beam 
bunches crossed at the IR. These, combined with asym- 
metry estimates from beam background dominated sam- 
ples, resulted in 0.003 uncertainty in All- 

The bias toward hard fragmentation processes caused 
by the HT trigger requirement was simulated, as were 
possible biases introduced by jet reconstruction and jet 
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FIG. 3: The longitudinal double-spin asymmetry All in 
p+p^ jet+X at ^fs = 200 GeV versus jet pt- The uncertain- 
ties on the data points are statistical. The gray band indicates 
the systematic uncertainty from the beam polarization mea- 
surement, and the hatched band the total systematic uncer- 
tainty. The curves show predictions based on deep-inelastic 
scattering parametrizations of gluon polarization @, [3- 



Pt resolution. The resulting pt dependent shifts in All 
were estimated with the polarized parton distributions of 
Ref. [2(J. Their total is estimated to be less than 0.009. 

Analyses with randomized proton beam helicity config- 
urations and other cross checks including parity violating 
single-spin asymmetries showed the expected statistical 
behavior, thus indicating no evidence for beam bunch to 
bunch or fill to fill systematics in All- 

The curves in Figure [3] show theoretical evaluations 0, 
[Lsj at p — px for the commonly used polarized parton 
distributions of Ref. J2(|. They shift by less than 0.003 
(0.017) at p T = 5.6 (15.7) GeV/c if fi is varied by a fac- 
tor of 2. The polarized parton distributions are based 
on a best fit to polarized inclusive DIS data, the so- 
called GRSV-standard gluon polarization distribution, 
and on assumptions of (i) a vanishing gluon polarization 
Ag(x, Qq) — 0, and (ii) maximally positive or negative 

gluon polarization, Ag(x, Qq) — ±g(x, Qn). at the initial 
">2 — O/i n^\r 2 Ir2 



scale Qq = 0.4 GeV /c in the analysis [201 ) . Alternative 
parametrizations [2l| are within the range spanned by 
(ii). Our data fall below the Ag(x, Qq) = g(x, Qq) eval- 
uation (x 2 /ndf ~ 3) and are consistent with the other 
evaluations, in qualitative agreement with Refs. [H, 0, [f|. 
The results thus disfavor large and positive gluon po- 
larization, proposed [22j originally to explain the small 
quark spin contribution to the proton spin. 

In summary, we report the first measurement of the 
longitudinal double-spin asymmetry All for inclusive 
jets with transverse jet momenta of 5 < pt < 17GeV/c 
produced at midrapidity in polarized proton collisions at 
y/s = 200 GeV. The jet cross section was determined 
for 5 < p T < 50 GeV/c and is described by NLO pQCD 
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evaluations over 7 orders of magnitude. The asymmetries 
All are consistent with NLO pQCD calculations utiliz- 
ing polarized quark and gluon distributions from inclu- 
sive DIS analyses, and disfavor at 98% C.L. large positive 
values of gluon polarization in the polarized nucleon. 
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